Kondo effects in a triangular triple quantum dot with lower symmetries 
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We study the low-energy properties and characteristic Kondo energy scale of a triangular triple 
quantum dot, connected to two non-interacting leads, in a wide parameter range of a gate voltage 
and distortions which lower the symmetry of an equilateral structure, using the numerical renormal- 
ization group approach. For large Coulomb interactions, the ground states with different characters 
can be classified according to the plateaus of = (<5 e — S )(2/ir), where S c and 8 are the phase shifts 
for the even and odd partial waves. At these plateaus of 0, both and the occupation number 
TVtot = (Sc + S )(2/n) take values close to integers, and thus the ground states can be characterized 
by these two integers. The Kondo effect with a local moment with total spin S — 1 due to a Nagaoka 
mechanism appears on the plateau, which can be identified by ~ 2.0 and AT to t ~ 4.0. For large 
distortions, however, the high-spin moment disappears through a singlet-triplet transition occurring 
within the four-electron region. It happens at a crossover to the adjacent plateaus for ~ 0.0 and 
~ 4.0, and the two-terminal conductance has a peak in the transient regions. For weak distor- 
tions, the SU(4) Kondo effect also takes place for iVtot — 3.0. It appears as a sharp conductance 
valley between the S — 1/2 Kondo ridges on both sides. We also find that the characteristic energy 
scale T* reflect these varieties of the Kondo effect. Particularly, T* is sensitive to the distribution 
of the charge and spin in the triangular triple dot. 

PACS numbers: 72.10.Fk, 73.63.Kv 



I. INTRODUCTION 

The triangle is the simplest polygon, and has a closed 
loop which plays an important role on various fascinating 
phenomena in the condensed matter physics. The closed 
path in a metal and semiconductor allows the electrons to 
move around the loop, and causes a quantum-mechanical 
interference effects, such as an Aharanov-Bohm (AB) 
effect The closed path consisting of the odd- number 
of links also causes frustration, which leads to resonating 
valence bonds for some anti-ferromagnetic systems^ 

Furthermore, the interplay between the strong correla- 
tion and the interference effects caused by the triangular 
structure has also been one of the topics of the current in- 
terests in different fields of the condensed matter physics. 
For instance, the single triangle is also a fundamental unit 
of the triangular and kagomc lattices. In these systems 
the geometrical frustration affects significantly the mag- 
netic properties and the behavior at the Mott-Hubbard 
metal-insulator transitions^ Another interesting exam- 
ple is the triangular trimer of Cr atoms placed upon an 
Au surface^— and this system is expected to show a non- 
Fermi-liquid behavior due to the multi-channel Kondo 
effects 

Recently, the triangular triple quantum dot (TTQD) 
has been experimentally realized and intensively stud- 
ied using various systems, such as AlGaAs/GaAs 
heterostructure*^— and self-assembled InAs quantum 
dotpi£ Theoretically, the TTQD has been shown to 
demonstrate various types of the Kondo effects^ 7 - - — as 



well as the AB effect The closed path makes the 
TTQD different from a linear quantum-dot chair>2£~— 
and the other three-level systems^ 3 - - — One of the most 
interesting points is that the appearance of a local mo- 
ment with total spin S = 1, at the filling where one ad- 
ditional electron is introduced into half- filling! 17 ' 39 This 
is caused by a Nagaoka ferromagnetic mechanism for the 
electrons moving around the triangular structured The 
S = 1 moment shows a Kondo behavior when the leads 
are coupled to the quantum dotspi 7 - - — Another interest- 
ing point is that the SU(4) Kondo effect takes place at 
half-fillingfi 7 - - — in the case where the ground state has a 
4-fold degeneracy caused by the orbital and spin degrees 
of freedoms. The TTQD has provided a new variety to 
the SU(4) Kondo effect, which had been studied for the 
double-dot systems4i~— 

The number of leads connected to the TTQD also af- 
fects significantly the Kondo behavior. This is because 
whether or not the local moment can be screened depends 
on the relation between the dimension of the Hilbcrt 
space for the local moment and the number of conducting 
channels^ The low-temperature properties of the TTQD 
have been studied, so far, by several theoretical groups, 
for the configurations with one^ two^ ' 20 ' 22 and three 
leadsd These studies complement each other the wide 
parameter space of the TTQD. Zitko et a&- and Mitchell 
et al^- studied the Kondo effect at half-filling in some sit- 
uations, but the dependence on the electron filling was 
not examined. Vernek et a&- examined the gate voltage 
td dependence in a wide range of the electron filling, but 
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the parameters used were confined to a region of small 
interaction U and a large dot-lead coupling T where the 
Kondo behavior is still rather smeared. 

We have studied the Kondo behavior of the TTQD 
away from half-filling in the series of the works Our 
research in the early staged started with a theoretical 
observation of a two-stage Kondo screening of the S = 1 
Nagaoka high-spin at four-electron filling and a sharp 
conductance dip caused by the SU(4) Kondo effect at 
half-filling (with three electrons), appearing in the gate- 
voltage dependence. The precise features of these Kondo 
effects have been clarified further in the previous paper ^ 
for the parameter values which cover the weak and strong 
couplings with respect to both U and T. We have also ex- 
amined the effects of the perturbations which break the 
full symmetry of the equilateral triangle,— as the real 
TTQD systems have some deviations from the regular 
structure in most of cases. Our results, obtained with 
the Wilson numerical renormalization group (NRG) ) 45 i 46 
have shown that the conductance dip typical of the SU(4) 
Kondo effect in the TTQD is quite sensitive to the per- 
turbation, while the 5 = 1 Kondo behavior is robust* 1 ^ 

The distortions of the triangular structure discussed in 
the previous paper, however, were still relatively small, 
so that the overall features of the Kondo effect in the 
TTQD have not yet been fully revealed and much re- 
mains to be explored. Particularly, a singlet-triplet tran- 
sition between a local singlet and the Nagaoka high-spin 
state occurs in the isolated TTQD cluster for large defor- 
mations, and this transition will affect the Kondo behav- 
ior at four-electron filling. Furthermore, the behavior of 
the conductance dip due to the SU(4) Kondo effect also 
needs to be clarified in more detail. 

The purpose of the present work is to provide a com- 
prehensive overview of the Kondo effect in the TTQD 
and to study the effects of large distortions. Specifically, 
we examine two different types of distortion: (i) an ir- 
regular inter-dot coupling, and (ii) an inhomogeneity in 
the level position of the quantum dots. We calculate the 
phase shifts, S c and S a , for the even and odd partial waves 
of the renormalized quasi-particles, in a wide parameter 
region of the gate voltage td and the distortions. These 
two phase shifts determine the ground state properties of 
the TTQD connected to two leads. 

In the parameter space for large U and small T, we 
find plateau with the integer values of the phase differ- 
ence O = (S c — S )(2/n), and at each plateau the oc- 
cupation number given by the Friedel sum rule N to t = 
(S c + <5 Q ) (2/ 7r) also approaches to an integer. These 
plateaus, therefore, can be classified with the two integer 
set (iVtot, ©), and each plateau corresponds to the ground 
state realized in each parameter region. For instance, 
the plateau for the 5=1 Kondo region can be labelled 
as (iVtot, 6) ~ (4.0, 2.0). The singlet-triplet transition 
emerges as a steep rise in O to the adjacent plateaus 
with the label (N tot , 6) ~ (4.0, 0.0) and (4.0, 4.0), situ- 
ated in the regions of a large distortion. Therefore, the 
two-terminal conductance g s shows a peak of the unitary- 
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FIG. 1. Triangular triple quantum dot in (a) series and (b) 
parallel configurations. The dot which has no direct connec- 
tion to the leads is referred to as the apex site, and labeled as 
i — 2 in the text. 



limit value g s = 2e 2 /h in the middle of the rise. We also 
find that the SU(4) Kondo behavior appears in the pa- 
rameter space along the contour line for O = 2.0, which 
traverses the middle of the steep rise in O between the 
plateaus with O ~ 1.0 and O ~ 3.0, for 7V to t — 3.0. 

We also estimate the characteristic energy scale T* of 
the Kondo screening in the wide parameter region, from 
the flow of the low- lying excitation energies in the NRG. 
The energy scale T* depends strongly on the local charge 
distribution in the TTQD. The screening is protracted 
significantly in the case where the partial component of 
the local moment becomes large at the apex site, which 
is located away from the leads as shown in Fig. Q] (a). 
This is because that the tunneling processes of the con- 
duction electrons from the leads to the apex site tend to 
be suppressed in the intermediate states on the other two 
sites. 

The paper is organized as follows. We describe the 
model and the formulation in Sec. [TTJ Some character- 
istics of the TTQD, seen already in the non-interacting 
case of U = 0, are summarized in Sec. IIIII Then, the 
molecular limit T — > for finite U is considered in Sec. 
IIVI in order to see the basic features of the local charge 
and spin states of the TTQD. The NRG results for the 
ground-state properties are shown in Sec. [V] The results 
for the characteristic energy scale T* are presented in 
Sec. ED A summary is given in Sec. IVHl 



II. FORMULATION 
A. TTQD connected to two non-interacting leads 

We consider a three-site Hubbard model on a triangu- 
lar cluster as a model for the TTQD. The cluster is con- 
nected to two non-interacting leads on the left (L) and 
right (R) as illustrated in Fig.[T](a). The Hamiltonian is 
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given by 

H = ^dot ^ckrt ^mix + Ulead , 

Nn 



(1) 
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^dot = X "'Mt "<W, "<Ma = d L d i*> ( 3 ) 



i = l 



H mix = X ( 4^ + v R d f ND , a C Ra + H.c. ) , (4) 



Klead= X X £ 



(5) 



v—L^R ha 



Here, d\ a creates an electron with spin a at the i-th dot, 
td,i the onsite potential, U the Coulomb interaction, and 
the number of the dots is given by Nd = 3. The hopping 
matrix element t%j between the dots is chosen to be real 
and positive (ty > 0). The dots labelled by i = 1 and 
z = 3 arc directly coupled, respectively, to the left and 
right leads via the tunnelling matrix element v L , R . The 

coupling causes the level broadening of F L / R = ^P v l/r> 
with p the density of states for the conduction band de- 
scribed by ek, and we will take T to be a constant as- 
suming a wide flat band. The conduction electrons are 



described by the operators ct and C. 



E 



k ^kh 



r/VN. 



In the present work, we consider the case that the sys- 
tem has an inversion symmetry choosing T L = T R (= T) , 



namely v l 



(= v), e dA = e di3 (= e d ), t 12 = t 23 (= t) 



and ti3 (= £'). We shall refer to the dot which has no 
direct connection to the leads as the apex site, and will 
use a notation ed,2 = fapcx for i = 2. We also choose the 
Fermi energy Ef as the origin of the energy Ep = 0. 



B. Phase shift, conductance and local charge 

Charge transfer between the dots and leads makes 
the low-energy states of the whole system a local Fermi 
liquid, which can be described by renormalized quasi- 
particles. Specifically, in the inversion symmetry case the 
two phase shifts, S c and d , for the even and odd partial 
waves become the essential parameters which character- 
ize the ground state [see also Appendix [A"] . At zero tem- 
perature T = 0, the series conductance g s for the current 
flowing through the two-channel configuration shown in 
Fig. H] (a), and the total number of electrons N tot in the 
dots can be expressed in terms of these phase shifts ; 19 ' 47 



g s = g sin (S e - 6 Q ) , 

Nn 



(6) 
(7) 



where g = 2e 2 /h. Note that the sum N to t and the dif- 
ference Q = (5 e — 5 )(2/n) between the two phase shifts 
link directly to the ground-state properties in the series 
configuration. Specifically, becomes a more natural 
measure for classifying the parameter space than <7 S for 
the quantum dots consisting of more than three local or- 
bitals Np> > 3. This is because the phase difference can 
be greater than n, for instance, it takes a value in the 
range < <5 C — S Q < 2tt for Nd = 3. 

The parallel conductance g p for the current flowing 
along the horizontal direction in the four-terminal geom- 
etry, shown in Fig. [T](b), can also be deduced from these 
two phase shifts 5 a and S Q defined with respect to the 
series configuration, 
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g (sin 2 5 e + sin 2 5 ) 



(8) 



The even and odd channels contribute to the parallel con- 
ductance separately with no cross terms which would 
represent interference effects. Note that in the case 
where the series conductance reaches the unitary-limit 
value g s = 2e 2 /h, namely at <5 e — 5 = (n + l/2)7r for 
n = 0, ±1, ±2, . . ., the parallel conductance also takes 
the same value g p = 2e 2 /h which is the half of its maxi- 
mum possible value 4e 2 jh. 

The phase shifts for the interacting case can be ex- 
pressed in terms of the renormalized hopping matrix el- 
ement tij for the quasi-particles ; 19 ' 48 
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(9) 



Here, S^(cj) is the self energy due to the Coulomb inter- 
action n^ ot , defined in Appendix EJ For the TTQD, the 
renormalized matrix elements expressed in the form 
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(10) 



The explicit form of the phase shifts can be obtained by 
solving the scattering problem of the renormalized quasi- 
particles, or equivalcntly from the Dyson equation given 
in given in (|A2[) . as 



cot dp = 



£rf — il — 2£ 2 /e ap ex j. £d + t' 



-, C0t5 o = 



r 



(ii) 



Specifically, the zero points of g s can be determined by 
the condition between the renormalized parameters 



F t' 4- f 2 — 



(12) 



which follows from the relation cot 5 e = cot 5 . Similarly, 
g s takes the unitary-limit value g s = 2e 2 /h in the case of 
cot 5 C = — l/cot<5 , which corresponds to the condition 

(e d - ? - 2?/? apcx ) (? d + ?) + T 2 = . (13) 



4 



III. EFFECTS OF DISTORTIONS IN THE 
NON-INTERACTING CASE 

We first of all discuss the level structure of an isolated 
cluster for T — in order to trace out the particular 
characteristics of the TTQD. We then calculate the con- 
ductances through the dots for U = to see how they 
reflect the level structure. These examples provide us 
with essential information for understanding the variety 
of forms which we will encounter in the wider parameter 
space. 



A. Level structure of the TTQD 

The one-particle energy levels for the non-interacting 
TTQD cluster which is described by H® ot are given by 



E. 



(l) _ e ap cx + £d — t 



e,± 



E^ = e d + t' 




2t 2 , 
(14) 
(15) 




-3-2-10 1 2 3 



FIG. 2. (Color online) Ground-state phase diagram of the 
isolated TTQD for Y -¥ and U = is plotted in (a) an e<i 
vs t' plane for e apex = e^, and in (b) an vs Ae plane for 



t' = t. Here, Ae = e a 



- £d. The phase boundaries are given 



by the contours for the energy levels, crossing the Fermi level: 
E { 1] = (solid line), E^l = (dashed line), and E^]_ = 
(dot-dashed line). The label (N tot , Q) is assigned for each 
region divided by these lines, and it represents the occupation 
number of the TTQD N tot , and 6 = (5 e - S )(2/w), which 
only for U = coincides with the difference in the occupation 
number of the even-parity levels and that of the odd-parity 
level. 



Here, E^ and E^ are, respectively, the energy for the 
cigenstatcs with the even and odd parities [see also Ap- 
pendix [Bj. Among the three eigcnstates, the one with 

the energy E^_ is the lowest for t > and t' > 0. The 

excited states become degenerate, E^ — E X \ for an 
equilateral triangle with e ape x = e<2 and t' = t. The de- 
generacy is lifted as the symmetry is lowered by a site 
diagonal distortion Ae = e apcx — e c i and also by an off- 
diagonal distortion At = t' — t. The first order correction 
is given by - E Q 1] ~ (2/3) Ae - (4/3) At. Thus, for 
Ae — 2 At > (< 0), the energy of the even excited state 
becomes larger (smaller) than that of the odd one. It 
reflects the fact that e ape x belongs to the even part of the 

basis, and the odd energy Eo increases with t' 

The number of the electrons which enter the TTQD is 
determined by the relative position of these levels with 
respect to the Fermi energy Ep (= 0). Figure [2] shows the 
phase diagrams of the ground state of the TTQD cluster 
for {7 = and r — s- 0. The boundaries are determined 
by the condition that the one-particle energy level crosses 



the Fermi energy: EP = (solid line), = (dashed 

line), and E^_ = (dot-dashed line). The left panel 
(a) is plotted as a function of td and t' keeping the site 
diagonal part uniform e apcx = e^. Similarly, the right 
panel (b) is plotted as a function of e c i and Ae keeping 
the inter-dot couplings uniform t' = t. Therefore, in the 
case of off-diagonal distortion, the cluster deforms from 
the regular triangle to a linear chain for < t' < t, and 
then for t' > t the coupling between the apex site and the 
other two becomes relatively weak as t' increases. The 
diagonal distortion e apcx ^ e^ affects directly the charge 
density in the apex site, and in Fig.[5](b) the contour for 



(i) 



the odd level becomes a vertical line because E^ does 
not depend on e apex . 

The occupation number A^totj varies discontinuously 
as an energy level crosses the Fermi energy. For fi- 
nite F, it can be deduced from the Friedel sum rule 
(S c + 5 )(2/n). As shown in Fig. [2l it takes the values 
Ntot = 0, 2, 4, and 6, depending on the region that is 
separated by the boundaries. The difference in the two 
phase shifts O = (<5 e — <5 )(2/7r) coincides, for U = 0, with 
Ncvcn — ^V ddj where N even and N oa \d are the occupation 
number for the even and odd levels, respectively. In the 
limit of r — > 0, it can take the values of O = 0, 2 and 4 
in the case of the TTQD, as shown in Fig. [2j For inter- 
acting electrons U ^ 0, however, there is no such general 
correspondence between Q and the charge difference in 
the even and odd subspaces, while the Friedel sum rule 
remains valid. This is because the Coulomb interaction 
T~0l ot breaks the charge conservation in each subspace, as 
seen explicitly in Eq. (|B10j) . 

We also see in Fig. [2] that the odd-parity level E^ 
(solid line) and the excited even-parity level E^]_ (dashed 
line) cross each other at t = t' and e apcx = ea, where the 
system has the full symmetry of the equilateral triangle. 
The crossing divides the region of four-electron occupa- 
tion into two different spin-singlet regions, which can be 
classified according to the values of 9. In the region with 
G = the highest occupied orbital is the odd-parity b\ 
orbital, while in the opposite side with = 4 the the 
even-parity orbital with energy E^ 1 } becomes the high- 
est occupied orbital. Note that the lowest even-parity 
orbital with energy E^l_ has already been occupied by 
two electrons in this area of the parameter space. 
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(a) e d /t (b) e d /t 

FIG. 3. (Color online) Series (a) and parallel (b) conduc- 
tances in the noninteracting case U = 0.0 as functions of ed/t 
and t'/t, for T/t = 0.25 and e apex = e d - The values of the 
conductances are scaled by go = 2e 2 /h. In (a) the solid line 
denotes the contour for S c — 8 a = n, and the dashed lines are 
the contours for the values of ir/2 and 3n/2. 




(a) e d /t (b) e d /t 



FIG. 4. (Color online) Series (a) and parallel (b) conductances 
in the noninteracting case U = 0.0 as functions of ed/t and 
Ae/t, for T/t = 0.25 and t' = t. Here, Ae = e apcx - e d . The 
values of the conductances are scaled by go = 2e 2 /h. In (a) 
the solid line denotes the contour for 8 C — S = n, and the 
dashed lines are the contours for the values of 7r/2 and 3n/2. 



B. Conductance for U — 

We next consider the noninteracting TTQD which are 
connected to the leads in a series or parallel configura- 
tions, as shown in Fig. [TJ In this case, the conduction 
electrons from the leads are scattered at the TTQD. The 
phase shifts 6 C and 6 , caused by the scattering, are given 
by Eq. (fTTj) , replacing the renormalized parameters there 
by the bare ones t, t' , and e ape x- The conductance can 
be deduced from these phase shifts through Eqs. ^ and 
([5]). or equivalently from the Green's function using Eqs. 
(fXjj) and (|A4]) given in Appendix [£] 

The series and parallel conductances in the non- 
interacting case U = are plotted in Fig. [3] as functions 
of ed/t and t'/t for T/t= 0.25 keeping the onsite poten- 
tial for the three dots to be the same Ae = 0. Similarly, 
in Fig.@]the conductances are plotted as functions of ed/t 
and Ae/t, taking the inter-dot hopping matrix elements 
to be uniform t' = t. The figures [3J and [4] can be com- 
pared, respectively, to the phase diagrams given in Fig. 
121(a) and (b). 

Both g s and g v are enhanced as the resonance levels 
which correspond to the one-particle energies defined in 
Eq. ([14"]) and (|T5l) cross the Fermi level. The series and 
parallel conductances show a similar behavior in most of 
the parameter regions. We can see, however, that they 
show a quite different behavior at the point t = t' and 
Ae = 0, where the series conductance vanishes while the 
parallel conductance takes the maximum possible value 
g p = Ae 2 /h for two conducting channels. At this point, 

the two one-particle levels -Eg 1 } and E^ cross the Fermi 
level simultaneously, and the phase shifts take the value 
5 C = 3ir/2 and S = ir/2. 

The solid line in Fig. [3] (a) and Fig. 2] (a), denotes 
the contour of the difference in the two phase shifts for 
the value S c — S Q = it. Thus, this contour corresponds 
to a zero line of the series conductance, and it means 
that destructive interference is most pronounced along 
this line. Similarly, the dashed lines in Fig. |3J (a) and 
Fig. @] (a) are the contours for S e — S Q = n/2 and 3ir/2, 



on which the two conductances show peaks of the same 
height g s = 2e 2 /h and g p = 2e 2 /h. 

We can also see in Figs. [3J and O that some conductance 
peaks are sharp and the others are relatively wide. Par- 
ticularly, the resonance peak for the excited even-parity 

level E^l, which corresponds to the dashed line in Fig. [2] 
(a) and (b), is much sharper than the other peaks. This 

is because the eigenstate for E^ 1 } has a large spectral 
weight at the apex site which has no direct couplings to 
the leads, and thus the hybridization with conduction 
band is suppressed. This feature can also be seen in the 
explicit expression for the spectral weight for the non- 
interacting TTQD is given in Eq. flB8j| in Appendix |U 
Conversely, the resonance width becomes large for the lo- 
cal states, the spectral weight of which is mainly on the 
other two dots coupled directly to the leads. 



IV. GROUND STATE IN THE MOLECULAR 
LIMIT: T = AND U^0 

The model can be solved also for finite Coulomb inter- 
action U in the molecular limit T — > i 19 i 39 In this case 
the TTQD is disconnected from the leads, and described 
by the Hamiltonian, 

"Hdot = ^dot + "^dot ■ (16) 

The eigenstates of "Hdot determine the high-energy prop- 
erties of the system, particularly the properties of the 
local excitations near the quantum dots. It gives us a 
knowledge how the parameter space could be classified; 
such a classification relates to the fixed points of the 
renormalization group. In this section we examine the 
effects of the distortions on the local electronic states in 
the interacting case. 

Figures [5] and [6] show the phase diagram of the ground 
state of the isolated TTQD for T — s- 0. The Coulomb 
interaction is chosen to be (a) U/(2irt) — 0.2, and (b) 
U/(2irt) = 1.0. These figures can be compared with 
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(a) e d /U (b) e d /U 

FIG. 5. (Color online) Ground-state phase diagram of the 
TTQD for T = is plotted in a e d /U vs t'/t plane for Ae = 0. 
The Coulomb interaction is chosen to be (a) U/(2nt) = 0.2, 
and (b) U/(2nt) = 1.0. The occupation number iVtot is shown 
for each region. The total spin S is 1/2 for odd N to t, and 
S — for even N to t, except it becomes S = 1 in the middle 
of three iVtot = 4 regions. Along the horizontal line in the 
iVtot = 3 region, the ground state has an SU(4) symmetry 
due to the parity and spin degeneracies. The eigenstate has 
an even (odd) parity above (below) this horizontal line. 




(a) e d /U (b) e d /U 

FIG. 6. (Color online) Ground-state phase diagram of the 
TTQD for T = is plotted in a e d /U vs Ae/U plane for t' = t. 
The Coulomb interaction is chosen to be (a) U/(2nt) — 0.2, 
and (b) U/(2wt) = 1.0. The occupation number Ntot is shown 
for each region. The total spin S is 1/2 for odd N to t, and 
S — for even JVtot, except it becomes S = 1 in the middle 
of three iVtot = 4 regions. Along the horizontal line in the 
ATtot = 3 region, the ground state has an SU(4) symmetry 
due to the parity and spin degeneracies. The eigenstate has 
an odd (even) parity above (below) this horizontal line. 



the phase diagrams in the non-interacting case shown 
in Fig. [2J We can see that the eigenstate with an odd- 
number of electrons (JVtot = 1, 3, and 5) and total spin 
5 = 1/2 becomes a ground state due to the Coulomb in- 
teraction. The odd-number electron regions emerge be- 
tween the even-number electron regions in the parameter 
space, and become wider as U increases. We can also see 
that three horizontal border lines appear in the phase 
diagrams for U > 0, and along each line a level crossing 
takes place between the two different eigcnstates with the 
same occupation number. 

The ground state for iVtot = 3 is separated by one of 
these horizontal lines at t' = t in Fig. [5j and similarly 
by the one at Ae = in Fig. [S] The ground state has 
an SU(4) symmetry along this border due to the orbital 
degeneracy caused by the symmetry of the equilateral tri- 
angle and the spin degeneracy. Away from this horizontal 
line, the distortions lower the equilateral symmetry, and 
lift the orbital degeneracy. An even-parity (odd-parity) 
state becomes the ground state for t' > t (t' < t) in 
Fig. [5] Correspondingly, in Fig. [5] the ground state is an 
even-parity (odd-parity) state for e apox < e d (e apox > e d ). 
Note that there are some similarities between the phase 
diagrams in Fig. [5] and Fig. [6) the features seen for t' > t 
(t' < t) are similar qualitatively (graphically) to those for 
Ae < (Ae > 0). This is because the two types of the 
distortion, t'/t and e apcx — e^, lift the degeneracy in an 
opposite way, as mentioned in the above with Eqs. (| 14[) 
and (fl5j) . 

The Coulomb interaction also causes the high spin S = 
1 ground state seen in Figs. [5] and [5] in the middle of 
the Ntot = 4 regions where the TTQD has one extra 
electron introduced into the half-filled cluster. The S = 
1 region evolves in the parameter space from the level 
crossing point for {7 = 0, seen in Fig. [5] at the point 
of t' = t and e ap ex = £d- The degeneracy at this level 



crossing point is lifted for infinitesimal {7, and the high- 
spin state evolves continuously, as U increases, to the 
Nagaoka ferromagnetic state which is usually defined in 
the large U limit. For large distortions, however, the 
transition to a singlet ground state takes place on the 
horizontal lines, running on the top and bottom of the 
S = 1 region in Figs. [5] and HI 

The isolated TTQD which is not connected to the leads 
has a local moment of S = 1/2 for odd-number fillings, 
and a high-spin S = 1 for Ntot = 4, as mentioned in 
the above. In the case where two leads are coupled to 
the cluster, however, the local moment is screened even- 
tually at low energies by the conduction electrons tun- 
neling from the leads, and the ground state of the whole 
system becomes a spin singlet. We show the results of 
the ground-state properties of the TTQD connected to 
the leads in the next section, and then discuss also the 
characteristic energy scale of the Kondo screening in Sec. 

ED 



V. NRG RESULTS FOR GROUND-STATE 
PROPERTIES OF THE INTERACTING TTQD 

We now consider an interacting TTQD coupled to two 
leads via tunneling matrix elements % m j X defined in Eq. 
(H|). In this case the phase shifts S c and S play an cen- 
tral role on the low-energy properties. The effects of 
the Coulomb interaction enter through these two phase 
shifts, which can be expressed in terms of the renormal- 
ized parameters for the quasi-particlcs of the local Fermi 
liquid, as described in Eq. ([11]) . 

We have calculated the many-body phase shifts using 
the NRG method) 19 ' 49 and have deduced the conductance 
and the occupation number of the TTQD at zero tem- 
perature from the phase shifts^ using Eqs. (jl])-©. In 
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our calculations, the ratio of the inter-dot hopping ma- 
trix element t and the half width of the conduction band 
D, defined in Appendix [C] is chosen to be t/D = 0.1. 
The iterative diagonalization has been carried out by us- 
ing the even-odd basis, described in Appendix For 
constructing the Hilbert space in each NRG step, instead 
of adding two orbitals from even and odd orbitals simul- 
taneously, we add one orbital from the even part first and 
retain 3600 low-energy states after carrying out the diag- 
onalization of the Hamiltonian. Then, we add the other 
orbital from the odd part, and again keep the lowest 3600 
eigenstates after the diagonalization. The discretization 
parameter is chosen to be A = 6.0, which has been con- 
firmed to reproduce the noninteracting results with a suf- 
ficient accuracy.— ~— In the following, we set the strength 
of the Coulomb interaction to be U/(2nt) = 1.0, which 
is adequate for observing typical results caused by U, as 
seen in Figs. [5] and |6] We have carried out some calcula- 
tions changing U and T for the equilateral triangle in the 
previous workJ^ Our results have clarified how the value 
of U affects the width of the Kondo ridges and the Kondo 
energy scale. Furthermore, a large Y smears the elec- 
tronic structure of the TTQD origin. Through these ob- 
servations, we have confirmed that the characteristic fea- 
ture of the Kondo behavior can be seen clearly for typical 
a parameter set of U/{2irt) = 1.0 and Y/t = 0.12, 0.25. 



A. Off-diagonal distortions: t' ^= t 

We discuss in this subsection the transport properties 
in the presence of the off-diagonal distortion t' ^ t keep- 
ing the site-diagonal potential uniform Ae = 0. The 
effects of the diagonal distortion e ape x ^ £d are examined 
in the next subsection IV Bl 



Total Charge in TTQD 




(b) 



-0.5 




FIG. 7. (Color online) NRG results for the occupation num- 
ber JVtot for U/(2irb) = 1.0, Y/t = 0.25, and Ae = 0. 
The left panel (a) shows JVtot as a function of Ed/U for 
t'/t = 0.0, 0.5, 1.0, . . . , and 3.5. The right panel (b) shows 
JV to t in the ed/U vs t'/t plane. The dotted lines in (b) are the 
contours for JVtot = 0.5, 1.0, 1.5, . . . , and 5.5 (in steps of 0.5 
from the right to the left). 



contours for JV tot = 0.5, 1.0, 1.5, . . . , and 5.5 (in steps of 
0.5 from the right to the left). Note that this figure can be 
compared with Fig.[5](b) where JVtot for r — ► is shown. 
We can see in Fig. 0(b) that the contours of JVtot for half 
integers (0.5, 1.5, . . . , and 5.5), which are shown with the 
thicker dotted lines, follow almost faithfully the phase 
boundaries between the different charge states for Y —> 
shown in Fig. [5] (b). The electron filling JVtot changes 
rapidly near these contours for the half integers. This can 
be seen explicitly in Fig. [7] (a), and it reflects the fact that 
r is much smaller than the inter dot matrix elements t 
and the Coulomb interaction U in the present case. From 
these observations, we see that the charge distribution in 
the plateaus regions is almost completely determined by 
the high energy states, and it can be approximated by the 
one in the limit of Y — > 0. The low-lying energy states 
are required, however, to describe correctly the transport 
properties and the conduction-electron screening of the 
local moment of the TTQD. 



1. Local charge: JVtot 



^(S e + S ) for t^t' 



Figure [7] shows the NRG results of the occupation num- 
ber JVtot for Y/t = 0.25, and U/(2id) = 1.0. In (a), the 
results are plotted as a function of ed/U for several of 
values of t'/t (= 0.0, 0.5, . . . , and 3.5, in steps of 0.5). 
We can see clearly that the plateaus emerge near inte- 
ger values of JV tot due to the Coulomb interaction, es- 
pecially the one for JVtot — 3.0 becomes almost flat for 
large t'/t > 3.0. These results show that the coupling 
strength Y = 0.25t is small enough to distinguish the dif- 
ferent charge states for U = 2nt. The plateau for the 
five-electron filling emerges due to the distortion and it 
becomes wider as t'/t deviates from 1.0. These features 
are consistent with that for the isolated TTQD discussed 
in the above. 

We have carried out the calculations more densely for 
a number of points in the parameter space than those 
presented in Fig. [7] (a), and the results are plotted in the 
ed/U vs t'/t plane in Fig. [3(b). The dotted lines are the 



2. Series and Parallel Conductances for t ^ t' 

The NRG results for the conductances at zero temper- 
ature are shown in the ed/U vs t'/t plane in Fig. [5] for 
U/(2nt) = 1.0, Y/t = 0.25, and Ae = 0. In order to 
show more clearly the overall features, we have provided 
two types of the plots seen from different points in the 
parameter space for each of the conductances. The series 
conductance g s is plotted in (a) and (c). Similarly, the 
parallel conductance g p is shown in (b) and (d) . For com- 
parison, the phase boundary for Y — > given in Fig.[5](b) 
is also superposed onto Fig. [8] (b) with the dashed lines. 
We see that the feature of the conductances reflects the 
occupation number in each of the regions in the param- 
eter space. Note that the ground state becomes a spin 
singlet in the whole region of the parameter space due to 
the screening by the conduction electron. 

We can also see in Figs. [5] (a) and (c) that typical 
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Series Conductance 



Parallel Conductance 




FIG. 8. (Color online) Series (a) and parallel (b) conduc- 
tances for U/(2%t) = 1.0, F/t = 0.25 and Ae = are plotted 
as functions of ed/U and t' /t. In (a), the solid line is a con- 
tour for the phase-shift difference 5 C — S = n, and the dotted 
lines are also the contours for 8 C — 8 = n/2 and 37r/2. In 
(b), the dashed lines denote the phase boundaries in the iso- 
lated limit r —¥ 0, corresponding to the ones shown in Fig. 
[S] (b) . Lower panels: (c) and (d) are the surface plots of the 
series and parallel conductances corresponding to (a) and (b), 
respectively. 



Kondo ridges for the series conductance with g s ~ 2e 2 /h 
emerge for odd-number fillings N to t — 1.0, 3.0, and 5.0. 
Furthermore, both g s and g p almost vanish for even- 
number fillings iVtot — 0.0, 2.0, 4.0 and 6.0 except for 
the 5 = 1 Kondo region. Particularly, the behavior at 
small fillings iVtot S 2-0, for ed/U > —1.0, can be ex- 
plained simply by the 5 = 1/2 Kondo effect due to the 

lowest single molecular orbital of E^,_. Therefore, the 
characteristic features of the TTQD appear in the region 
°f e d ^ —0.2U, where the two excited levels Eo and 
-Ee + are partially filled. 

The solid line in Fig. [8] (a) denotes the contour for the 
difference in the two phase shifts corresponding to the 
value S c — S = n. Along this line, the series conductance 
becomes exactly zero. Specifically, in the three-electron 
region, for —0.8 < ed/U < —0.3, this contour runs near 
the horizontal line for t'/t = 1.0 where the triangle has 
the equilateral symmetry. The contour line tilts slightly 
from the horizontal line because the coupling to the two 
leads breaks the equilateral symmetry already at t' = t. 
This contour for 5 C — S Q = n appears in Fig.[5](c) as a very 
sharp valley of the series conductance. The Kondo ridges 
on each side of this valley have a different parity. Just 
at the bottom of the valley, the low-lying quasi-particle 
states for the even and odd channels become degener- 
ate, and the low-energy properties can be described by 
the SU(4) Fermi-liquid theoryJ^ Furthermore, along this 
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FIG. 9. (Color online) Ground-state properties at = 
— 0.6U, (a) conductances and (b) phase shifts (5 C ± S )(2/tt), 
are plotted in a narrow region of t' near t'/t = 1.0 for 
U/(2nt) = 1.0, F/t = 0.25, and ea. pex = e<i. In this parameter 
region, the occupation number is almost constant iVtot — 3.0. 



valley the two phase shifts are almost constant with the 
values, 5 ~ 7r/4 and <5 C ~ 57r/4, since the Coulomb in- 
teraction keeps the sum of the two to be 5 C + S ~ 37r/2 
in the three-electron region through the Friedel sum rule. 
Therefore, the parallel conductance does not change so 
much near this valley of the series conductance, keeping 
the value of g p ~ 2e 2 /h. 

In order to see the sharp SU(4) Kondo valley in more 
detail, the conductances and the phase shifts at = 
— Q.6U are plotted in Fig. [!|]as functions of t'/t. Partic- 
ularly, the two lines in Fig. [§] (a) correspond to a cross 
section of the surface plots given in Fig. [5] (c) and (d) in 
the middle of the three-electron region along the vertical 
direction. We can see in Fig. [S] (b) that the phase-shift 
difference S c — S increases with t'/t showing a kink, the 
value of which varies from 7r/2 to 37r/2 as t'/t increases, 
and taking the value of ir at t'/t ~ 1.0 in the middle of 
the transient region. This kink determines the structure 
of the series conductance valley seen in Fig. |H] (a). There- 
fore, the slope of the phase difference S c — 6 in the middle 
of the kink determines the width of the valley. Note that 
it is quite general to the local Fermi-liquid systems that 
the derivative of the phase shift with respect to the pa- 
rameters, such as the frequency and the external fields, 
plays an important role on the rcnormalization of some 
correlation functions. 

The 5=1 Kondo behavior can be seen for the four- 
electron filling in the diamond-shape region in Figs. [8] (a) 
and (b). The series conductance almost vanishes g s ~ 0.0 
in this region, while the parallel conductance is enhanced 
g p ~ 4e 2 /h despite an even-number electron filling. This 
contrast between g s and g p can be seen clearly, partic- 
ularly in Figs. [5] (c) and (d). We can also see in Fig. [5] 
(a) that the contour for S c — 5 = n is winding in the 
center of the diamond region near ed/U ~ —1.0. Such 
a bend is not seen in the noninteracting case, for which 
the contour varies monotonically as shown in Fig. [3] (a) . 
The contour lines of the phase shifts evolve, however, 
continuously from the non-interacting form. This is be- 
cause the ground state of the whole system evolves adi- 
abatically from a singlet described by a single Slater de- 
terminant to a correlated singlet described by the local 
Fermi-liquid theory for finite T. Note that the 5=1 mo- 
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FIG. 10. (Color online) Ground-state properties at e<; = 
— 1.043(7, (a) conductances and (b) phase shifts (5 C ±<5 ) (2/ 7r) , 
are plotted as functions of t'/t for U/(2nt) = 1.0, Y/t = 0.12, 



and e a 



e<j. In this parameter region, the occupation num- 



ber is almost constant iVtot — 4.0. Note that in the limit of 
r — > the Nagaoka state is the ground state for the isolated 
cluster for 0.30 < t'/t < 2.88. 



merit is screened at low temperature by the conduction 
electrons from the two leads via a two-stage screening 

1 7-1 9 

processes — 

The dotted lines in Fig. [5] (a) express the contours for 
S c — S = tt/2 (below the solid line) and 3tt/2 (above the 
solid line), on which the series conductance reaches the 
unitary-limit value g s = 2e 2 /h. One of the dotted lines 
on the right, at ed/U > 0.1, follows simply the Kondo 
ridge caused by the lowest orbital E^_. The other two 
lines pass on the top and bottom of the diamond of the 
S = 1 Kondo region. These two contours can be com- 
pared to the phase boundaries for the singlet-triplet tran- 
sition, seen in a narrow range of at t'/t ~ 0.30 and 
2.88 in Fig. [5] (b). In order to clarify the precise fea- 
ture of the corresponding crossover between the S = 1 
Kondo and non-Kondo singlet states, the conductance 
and phase shifts are plotted in Fig. [TU] as functions of 
t'/t, choosing the level position to be in the middle 
of the four-electron region at e<j ~ — 1.04C7. At each 
of the crossover points, near t'/t ~ 0.3 and 3.0, the se- 
ries conductance has a peak. The feature of these con- 
ductance peaks reflects the kink in the phase difference 
S c — S Q , the value of which varies from 0.2tt to I.Ott near 
t'/t ~ 0.3, and from I.Ott to 3.0tt near t'/t ~ 3.0. There- 
fore, the slope of these kink determines the width of the 
conductance of peak. Furthermore at the crossover re- 
gion, the electron occupation fluctuates slightly from 4.0 
as iVtot — 4.0 ~ ±0.5. In the 5 = 1 Kondo-singlet region 
situating between the two peaks of g s , the phase shifts 
are almost locked at S c ~ 3tt/2 and 8 a ~ tt/2, and thus 
the parallel conductance takes the value g p = 4e 2 //i. In 
one of the non-Kondo regions for t'/t > 3.2, the phase 
shifts approach to S c ~ 2tt and S — 0. The phase shifts 
take the values of S c ~ 1.2tt, and <5 — tt in the limit 
of t' — > 0, in the other non-Kondo region for t'/t < 0.2. 
Note that for Fig. QUI the coupling between the TTQD 
and the leads has been chosen to be T/t = 0.12, which is 
smaller than that (T/t = 0.25) for the previous figures, 
in order to see clearly the typical features of the narrow 
crossover regions. 




J.O -0.5 0.0 



FIG. 11. (Color online) Difference between the even and odd 
phase shifts @ = (5 e - S )(2/n) for U/(2nt) = 1.0, Y/t = 0.25 
and Ae — 0. The left panel (a) shows O as a function of td/U 
for the values of t'/t = 0.0, 0.25, 0.5, . . . , and 3.5 (in steps of 
0.25 from the bottom to the top). The right panel (b) shows 
G in the ea/U vs t'/t plane. The dotted lines in (b) are the 
contours for (S e — 5 )(2/ir) = 0.5, 1.0, 1.5, and 3.5 (in 
steps of 0.5 from the bottom to the top). 



3. Phase-shift difference: Q = —(8 C 



S ) for t ± t' 



The difference between the two phase shifts S c — S Q is a 
fundamental parameter that contains the essential infor- 
mation of the interference effects between the even and 
odd conducting channels. It affects the series conduc- 
tance, while each channel contributes independently to 
the parallel conductance, through the expressions given 
in eqs. © and ©. Specifically, peaks and dips of the se- 
ries conductance correspond directly to the kinks of the 
phase-shift difference <5 — S Q , as seen in Figs. [9] and [TU] It 
is also much easier for a numerical purpose to trace the 
kink structure of S e — S Q than to find directly the dips 
and peaks of g s . 

We also provide the NRG results for = (5 C — 5 q ){2/tt) 
in Fig. [TT]in order to clarify its behavior in the wide pa- 
rameter space. In the left panel (a), O is plotted as a 
function of ed/U for the values of t'/t varying from 0.0 
to 3.5 in steps of 0.25. Furthermore, Fig. [TT] (b) shows 
the results obtained in the ed/U vs t'/t plane: the dotted 
lines are the contours for the values of O varying from 0.5 
(bottom and right) to 3.5 (top) in steps of 0.5. We can 
see that there are several plateaus, or shelves, in these 
figures near the integer values of 9 ~ 1.0, 2.0, and 3.0, 
on which g s becomes almost transparent or zero. Fur- 
thermore, the occupation number also approaches to an 
integer value on each of these plateaus, and thus they 
can be classified according to a set of the two integers 
(iVtot, 6). For instance, in Fig. [TT](a), we can see a wide 
plateau which can labelled by (N tot , 0) ~ (3.0, 3.0) for 
-1-0 < e d /U < -0.3 and t'/t > 1.0. The height of the 
plateau, however, is still somewhat smaller than the ex- 
act integer 3.0. Such a deviation of the plateau height 
from an integer value decreases as T decreases. This has 
been confirmed explicitly for the equilateral triangle in 
the previous work [see Fig. 6 of Ref. OH . Furthermore, 
we can see another example for smaller T in the next 
section [see Fig. [TT]. 



10 



The feature of G in the parameter space can also be 
compared to the phase diagram for the isolated TTQD. 
Particularly, the contours for 9 = 0.5, 1.5, 2.5, and 3.5, 
which are shown with the thicker dotted lines in Fig. [11] 
(b), divide the parameter space in a similar way that the 
phase boundaries did in Fig. [5] (b). The contour lines, 
however, do not cross each other while the border lines 
for r = are crossing at some points. We can see in 
Fig. [TT] that the SU(4) Kondo effect is manifest in the 
parameter space as a sheer cliff at —0.8 < ed/U < —0.3 
near t' jt ~ 1.0. It also corresponds to the kink that we 
have seen in Fig. [5] (b). Between the bottom and top of 
the cliff the value of 8 varies from 1.0 to 2.9, respectively. 
The slope of the cliff determines the width of the SU(4) 
valley which corresponds to the contour line for = 2.0, 
running in the middle of the cliff. The Kondo ridges of g s 
on both sides of the valley can be classified according to 
the plateau value of <5 C — 6 ~ ir/2 or 37r/2, as the phase 
difference varies by 7r across the valley. 

The S = 1 Kondo region can also be seen as a diamond- 
shape plateau in Fig. [TT] (b), appearing at ~ —1.017 
and 0.3 < t' /t < 2.9. This plateau is characterized by 
the two parameters, ~ 2.0 and iVtot — 4.0. Thus the 
phase shifts are almost fixed at the value of S c ~ 37r/2 
and S ~ 7r/2 in this region. 



B. Diagonal distortions: e ap ox 7^ £d 

We next examine the effects of the diagonal distortion 
fapcx 7^ £d; keeping the intcr-dot hopping matrix elements 
uniform t' = t and taking the Coulomb interaction to 
be U = 2nt. In this subsection we choose the coupling 
between the leads and the TTQD such that T/t = 0.12, 
which is approximately a half of the one used for Figs. 
151 and [TT1 in the off-diagonal case. 




FIG. 12. (Color online) NRG results of the occupation num- 
ber JVtot for U/(2wt) = 1.0, for T/t = 0.12 and t' = t. The 
left panel (a) shows AT tot as a function of ed/U for values of 
for the values of Ae/U = -1.0, -0.9, -0.8, . . . , and 1.4 (in 
steps of 0.1 from the top to the bottom). The right panel (b) 
shows iVtot in the ed/U vs Ae/U plane. The dotted lines in 
(b) are the contours for N to t = 0.5, 1.0, 1.5, . . . , and 5.5 (in 
steps of 0.5 from the right to the left). 

to the left). Particularly, the thick dotted lines are the 
contours for the half-integer values; 0.5, 1.5, 2.5, 3.5, 4.5, 
and 5.5. We can see that these thick dotted lines almost 
follow the phase boundaries between the different charge 
states in the isolated TTQD, shown in Fig. [6] (b). The 
local charge 7V to t changes rapidly near these thick con- 
tours, and has a plateau of an integer value between the 
thick dotted lines, as can be seen explicitly in Fig. [J2] 
(a). Therefore the charge in the plateau regions is de- 
termined at high energy scale, and the sum of the phase 
shifts (2/7r)((5 e + S ) in the plateaus can be approximated 
reasonably by the value of iVtot in the T — > limit. How- 
ever, the transport properties at zero temperature are 
determined by each of the two phase shifts or the differ- 
ence between them, which are determined essentially by 
the low-lying energy states of the whole system including 
the leads. 



1. Local charge: N to t = — (S e + S ) for Ae 7^ 

Figure [T2l (a) shows the occupation number N tot in the 
TTQD as a function of ed/U for the values of Ae/U = 
-1.0, -0.9, -0.8, . . . , and 1.4 (in steps of 0.1 from the 
top to the bottom) . We can see clearly that the plateaus 
appear near integer values of iV to t- In the present case 
the coupling strength T is much smaller than U and t, 
so that the different charge states can be distinguished 
clearly. In other words, the crossover between two ad- 
jacent charge states becomes sharp, and thus the border 
can be determined reasonably by the middle point where 
.ZVtot takes a half-integer value. 

We have also carried out the calculations for a num- 
ber of parameter sets, more than the ones which are 
shown explicitly in Fig. [12] (a), and have plotted the re- 
sults in Fig. [T2] (b) in the e,i/U vs Ae/U plane. In this 
figure the dotted lines denote the contours for iV to t = 
0.5, 1.0, 1.5, . . . , and 5.5 (in steps of 0.5 from the right 



2. Series and Parallel Conductances for Ae 7^ 

The scries (a) and parallel (b) conductances are plotted 
in Fig.rT5]in the parameter space of ed/U and Ae/U . For 
comparison, the phase diagram for T —> given in Fig. 
[5] (b) is superposed onto Fig. [13] (b) with the dashed 
lines. We can see that the behavior of the conductances 
in this parameter space also reflects the feature of the 
phase diagram for the isolated TTQD. In the regions of 
the odd-number electron filling the both conductances 
g s and g p have the Kondo plateaus with the height of 
2e 2 /h. Furthermore, the S = 1 Kondo effect takes place 
in a trapezoidal region near Ae ~ 0.0 and ed/U ~ —1.0. 
In this region, the series and parallel conductances show 
a clear contrast, namely g s ~ while g p ~ 4e 2 /h. This 
feature is the same as what is observed in the case of the 
off-diagonal distortions. 

The solid line in Fig. [T3] (a) denotes the contour for 
5 C — S a = tt, on which the series conductance becomes 
zero. This contour runs across the region of the three- 
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(a) e„/U (b) e d /U 



FIG. 13. (Color online) Series (a) and parallel (b) conduc- 
tances for U/(2irt) = 1.0, T/t = 0.12 and t = t' are plotted as 
functions of ed/U and Ae/U. In (a), the solid line is a contour 
for the phase-shift difference S c — 5 a — tt, and the dotted lines 
are also the contours for the values of tt/2 and 3tt/2. In (b), 
the dashed lines denote the phase boundaries for the isolated 
TTQD, corresponding to Fig.[6](b). 

electron occupancy almost horizontally in an area with 
weak distortions Ae ~ 0.0. It associated with a sharp 
valley of the series conductance, which is typical of the 
SU(4) Kondo effect in the TTQD and is seen also for the 
off-diagonal distortions. The SU(4) symmetry is caused 
by the channel degeneracy restored along the line at low 
energies^ and the phase shifts take the values of 5 C ~ 
57r/4 and 5 Q ~ tt/A. We can also see that the contour for 
S c — S a = tt is deformed significantly, in the trapezoidal 
S = 1 Kondo region, from the non-interacting form which 
is a simple straight line shown in Fig. [4] (a). This could 
happen, however, continuously with increasing [/, as the 
ground state evolves adiabatically in the case that the 
quantum dots are coupled to the leads. 

The dotted lines in Fig. [T3] (a) denote the contours 
for 5 C — S = tt/2 (above the solid line) and 3tt/2 (be- 
low the solid line), on which the series conductance takes 
the unitary-limit value g s = 2e 2 /h. It should be noted 
that a long and very sharp ridge emerges in Fig. [13] 
(a) for the series conductance at Ae/U ~ —0.4 and 
—0.9 < ed/U < —0.6. This sharp ridge runs along the 
lower end of the trapezoidal S = 1 Kondo region, and 
reflects the singlet-triplet transition taking place in the 
isolated TTQD cluster for r — s- 0. 



2 

3. Phase-shift difference: O = — (5 C — S ) for Ae 7^ 

Figure [H] (a) shows the results of the phase-shift dif- 
ference O = (S e — 8 ){2/ix) as a function of ed/U for the 
values of Ae/U varying from —1.0 to 1.4 in steps of 0.1. 
We can clearly see that there are a number of plateaus 
near integer values of 6 ~ 1.0, 2.0, and 3.0. Specifi- 
cally, the height of the plateaus approaches very close to 
exact integers in the present case because the coupling 
between the leads and the TTQD is small. Although we 
can recognize that some of them, for instance, the ones 
near ~ 3.0, still deviate from an exact integer, these 
deviations can be controlled by tuning T to be smallJ^ 




FIG. 14. (Color online) Difference between the even and odd 
phase shifts 9 = (S e - <5o)(2/tt) for U/(2nt) = 1.0, T/t = 0.12 
and t' — t. The left panel (a) shows 6 as a function of ed/U 
for the values of Ae/U = -1.0, -0.9, -0.8, . . . , and 1.4 (in 
steps of 0.1 from the top to the bottom). The right panel (b) 
shows in the ed/U vs Ae/U plane. The dotted lines in (b) 
are the contours for (S c — 5 )(2/ir) — 0.5, 1.0, 1.5, . . . , and 
3.5 (in steps of 0.5 from the top to the bottom). 



In order to see the behavior of (5 C — 5 )(2/it) in the 
parameter space, the results arc plotted also in the ed/U 
vs Ae/U plane in Fig. Q3] (b). In this figure the dotted 
lines denote the contours for 0, particularly the thicker 
ones are the contours for the half-integer values: — 
0.5, 1.5, 2.5, and 3.5. Each of these thick dotted lines 
runs very closed to the phase boundaries for the isolated 
TTQD shown in Fig. [6] (b). These thick contours, as a 
whole, cover almost all the boundaries. These contour 
lines of 0, however, evolve continuously from the non- 
interacting forms as U increases. This is because the 
ground state of the whole system evolve adiabatically 
from a U — spin singlet to a correlated singlet described 
the local Fermi-liquid theory for finite T. 

The sharp SU(4) Kondo valley of the series conduc- 
tance corresponds to a steep cliff standing at —0.8 < 
td/U < -0.3 for small distortions Ae/U ~ 0.0 in Fig.[H 
(b). At the edge of this cliff, (6 C — 5 )(2/ir) varies by an 
amount 2.0 approximately. It varies from ~ 1.0 (for 
Ae > 0.0) to ~ 2.9 (for Ae < 0.0), taking the value 
of 5 C — 5 a = tt which corresponds to the zero point of g s 
in the middle of the cliff. The slope of this cliff deter- 
mines the width of the SU(4) valley, as that in the case 
of off-diagonal distortions. 

We can see another sharp cliff in Fig. [14] (b) at the 
bottom of the S = 1 Kondo region for Ae/U ~ —0.4 and 
— 0.9 < ed/U < — 0.6, where the singlet-triplet transition 
takes place for the isolated cluster. Between the top and 
bottom of this cliff, the phase difference varies rapidly 
from 2.0 to 4.0, which causes the sharp Kondo ridge of 
the series conductance, seen in Fig.Q]|](a). Note that the 
slope of the cliff determines the width of the sharp peak 
of g s . We can also see a narrow cliff due to the singlet- 
triplet transition, at the top of the S — 1 Kondo region 
for Ae/U ~ 0.8 and e d /U ~ -1.2. At this cliff, the value 
of changes from 0.0 to 2.0 approximately. 
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Characteristic Temperature 
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FIG. 15. (Color online) The characteristic energy scale T* 
in the e d /U vs t'/t plane for T/t = 0.25, U/(2irt) = 1.0 and 
Ae = 0. The results in the range of -20 < log 10 (T*/T) < 
are painted in the colors varying from blue to red. T* 
becomes very small in the purple region, and the solid lines 
there denote the contours for log 10 (T*/T) = —20, —40, and 
—60. The dashed lines denote the phase boundary in the limit 
of T — > 0, corresponding to Fig. [5] (b). 



VI. CHARACTERISTIC ENERGY SCALE 

The ground-state properties, discussed in the above, 
are determined by the behavior of the phase shifts for the 
quasi-particles. At finite temperature T, for instance, the 
structure of the plateaus and dips will be smeared grad- 
ually as T increases. Nevertheless, the corrections due to 
finite T can be still determined by the local Fermi-liquid 
theory for T < T*, namely at temperatures lower than 
the characteristic energy scale T*. Specifically, in the 
case the quantum dots have a local moment, T* can be 
regarded as the Kondo temperature such that the mo- 
ment is screened at T <C T* . The value of T* , however, 
depends sensitively on the parameters at each point of 
the parameter space. Therefore, the actual temperature, 
at which the Fermi-liquid behavior can be observed, is 
different depending on the region in the parameter space. 

In the NRG approach, the crossover from the high en- 
ergy to the low energy Fermi-liquid regime can be seen in 
the trajectory of the low-energy levels of the discretized 
Hamiltonian Hjy defined in Appendix [C] The trajectory 
evolves as the number of orbitals N of the conduction 
band increases. In the present work we have estimated 
T* through N* that is a particular value of N, at which 
the trajectory almost enters the low-energy regime, as 



(17) 



Therefore, T* is the typical energy scale of the low-lying 
excitations of a finite NRG chain with N = N* . The val- 
ues of T* determined in this way have some ambiguities, 
especially in the case where the crossover is gradual, and 
our definition tends to give a smaller value for the charac- 
teristic energy scale. Nevertheless, the relative variations 
of T* in the parameter space can be extracted reasonably 
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FIG. 16. (Color online) Charge distribution in the isolated 
TTQD, N a .i = X]<r( rl <MO') is plotted as a function of t'/t for 
(a) iV a o the apex site, and for (b) the even a± orbital described 
in Appendix E The average is taken with respect to some of 
the eigenstates of Udot = "Hdot + %dot> an d labelled by the 
occupation number, spin, and parity ("N3S05:even" denotes, 
"Ntot =3", S = 1/2, and an even parity). The parameters 
are chosen to be F = 0, U/(2nt) = 1.0, and e ap cx = £d- Note 
that the occupation number takes the form iVtot = N a ,o + 
N a ,i+N b ,x. 



as shown in Figs. fT5l and [l~8l We have also confirmed that 
T* determined in this way is consistent with the ones 
we had estimated from the entropy in the previous work 
for the equilateral triangle.— We will see below that the 
variations of T* can be understood roughly through the 
distribution of the charge and spin in the even and odd 
parity orbitals, described in Appendix [Bl 



A. T* vs off-diagonal distortions (t 1 7^ t) 

The results for T* are plotted in Fig. [T5] using a 
logarithmic scale, function of ed/U and t'/t for 

T/t = 0.25 and U/{2nt) = 1.0. The phase diagram of 
the isolated TTQD given in Fig. [5] (b) is also superposed 
onto Fig. Qj}] with the dashed lines, as a guide for the eye. 
We can see that the variation of T* in the parameter 
space also relates to the plateaus of [see Fig. [11] in Ap- 
pendix [B]. The energy scale T* is large in the case that 
the quantum dots have no local moment, and it becomes 
smaller when the TTQD has a local moment. 

We can see, nevertheless, T* is still relatively high in 
the regions of N to t — 1-0 and 5.0, in which the S = 
1/2 Kondo effect takes place. The S =1/2 moment for 
iVtot — 1.0 is caused by a single electron, which enters 
an even-parity state and stays mainly at the a\ orbital 
adjacent to the leads [see Fig. [21]. Thus the screening 
can be completed relatively easily in this case. Figure 
1161 shows the occupation number 7V a o (^V a i) m the ao 
(ai) orbital, as a function of t'/t, for the limit of V —> 0: 
note that ao corresponds to the apex site. The dashed 
line which is labelled "NlS05:even" shows the average 
with respect to the lowest eigenvector in the subspace 
with Ntot = lj S = 1/2, and an even parity. We can 
see that 7V Q i approaches to 1.0 as t'/t increases. In the 
opposite case, at t'/t = 0.0, the occupation number for 
ao and that for ai coincide N a Q = N a ± = 0.5, but still the 
charge and spin fluctuations are not suppressed because 
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FIG. 17. (Color online) Spectral functions in the noninter- 
acting case for (a) t'/t = 0.6 and (b) t' /t = 1.5. The solid 
line is the even component — Im [G a o{oj) + Gai(w)]/n, and the 
dotted line is the odd component — ImGbo{w)/n, defined in 
Appendix|B] The other parameters are chosen to be Ae = 0.0, 
e d = 0.0, and T/t = 0.25. 



these orbitals arc still at quarter filling. 

In the five-electron region for t'/t > 1.0, the eigen- 
vector |$^ d ) for the isola ted TTQD can be expressed 
in the form of Eq. (|B14[) . The local moment in this 
case stays at the b\ orbital, which is also adjacent to 
one of the leads [see Fig. [3T], and the conduction elec- 
trons can screen the moment through the usual kinetic 
exchange mechanism. For another five-electron region at 
< t'/t < 1.0, the eigenvector l^ivon) can be written in 
the form of Eq. (|B13|) . The averages iV a o and JV a i with 
respect to this state coincide with those with respect to 
the Nagaoka state l^^d) defined in Eq. (|B12[) , and the 
results are plotted in Fig. [Jj)] with the dot-dash line la- 
belled "N4Sl:odd". We can see that for < t'/t < 1.0 
a single hole with a spin 1/2 enters both of the even or- 
bitals, and for t'/t — >■ these orbitals approach to quarter 
filling in the hole picture. Therefore, the screening is not 
suppressed so much also in this five-electron region. 

The screening temperature becomes small in the three 
and four electron regions, namely —1.1 < < —0.3 in 
Fig. [HI We can see in the three electron region, however, 
T* is still much higher for t'/t < 1.0 than for t'/t > 1.0 
despite the local moment in the TTQD is S = 1/2 in both 
of the cases. For t'/t < 1.0, the ground state is an odd- 
parity state, and one of the three electrons enters the b\ 
orbital, and the other two electrons enter almost equally 
to the ao and a\ orbitals. This can be confirmed through 
the dotted line labelled "N3S05:odd" in Fig.QJU We can 
also see in Fig. [T7] (a) that the third electron enters the 
odd bi orbital for t'/t < 1.0 in the noninteracting limit. 
In contrast, for t'/t > 1.0 the ground state is an even 
parity state, and the solid line labelled "N3S05:even" de- 
notes average with respect to this state. We see that 
N a i approaches to 2.0 as t'/t increases, while the occu- 
pation of the apex site is almost unchanged N a o ~ 1.0, 
and thus the occupation of the b\ orbital is decreasing in 
this case. Therefore, the local moment is mainly due to 
the electron staying at the apex site. Thus the screening 
is protracted significantly because the charge and spin 
fluctuations are suppressed at the nearly filled a\ orbital, 
over which the conduction electrons come to screen the 
moment. We have also confirmed that along the sharp 



conductance valley caused by the SU(4) Kondo effect, at 
t'/t ~ 1.0 and -0.8 < e d /U < -0.3, the energy scale 
T* is enhanced due to the orbital degeneracy. Note that 
the variations in the spin and charge configurations near 
the SU(4) symmetric point becomes wider in the TTQD 
than the double dot. 

The properties of the local moment in the three- 
electron region also reflect the feature of the one-particle 
state which emerges as a peak of the conductances shown 
in Fig. [3] and also the corresponding spectral function is 
shown in Fig. [T7J Specifically for t'/t > 1.0 the orbital 
degeneracy is lifted such that E { 1] > E { c 1] + . Therefore, 
after two electrons occupy the lowest even-parity orbital 
with the energy E^_, the third electron enters the reso- 
nance state corresponding to the excited even-parity or- 
bital which appears in Fig. Q~7] (b) as the central peak. 
This state has a dominant spectral weight in the apex 
site, and the resonance width is narrower than T already 
in the noninteracting case [see also Eq. (|B8|) in Appendix 
IB] . It should also be noted that the width of this reso- 
nance determines T* in the noninteracting limit. For 
finite U, this peak may evolve into a Kondo resonance 
whose width is reduced further by the Coulomb interac- 
tion to the value of the order of T*. 

The wavefunction for the Nagaoka state has an 
odd parity, and the charge distribution of this state 
is shown in Fig. [16J with the dot-dash line labelled 
"N4Sl:odd". We can see that iV a0 and N al for the Na- 
gaoka state are similar to those for the three electrons 
state "N3S05:even" . The fraction of the local moment 
stays at the apex site, and this explains the reason why 
T* is small also in this case. One extra electron enters 
mainly the b\ orbital, and it provides half of the S = 1 
moment which can be screened at high temperature at 
the first stage of the two-stage Kondo screening^ 



B. T* vs diagonal distortions (e apcx 7^ e<i) 

The characteristic energy scale T* for the TTQD with 
the diagonal distortions is plotted in Fig. [15] using a 
logarithmic scale, as a function of ed/U and Ae/U for 
U/{2nt) = 1.0. The phase diagram of the isolated TTQD 
given in Fig. [SJ (b) is also superposed onto Fig. [TS] Note 
that the coupling between the leads and quantum dots is 
chosen to be T/t = 0.12, which is smaller than that we 
chose for the off-diagonal case. Therefore, the absolute 
values of T* become smaller than those in Fig. [J_5] We 
saw in the above that the Kondo screening is sensitive to 
the electronic structure of the TTQD for the off-diagonal 
distortions t' 7^ t. In order to see the results in a similar 
way, the charge distribution in the even and odd orbitals 
in the isolated TTQD for T — > is also plotted in Fig. 
[19] as a function Ae (= e apcx — e^). 

We can see that T* is suppressed also in the region of 
iVtot — 1-0 at Ae < 0.0, which corresponds to the area of 
ed/U > 0.3 at the right bottom of Fig. [15] In this param- 
eter region, the potential profile of the onsitc energy is 
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FIG. 18. (Color online) The characteristic energy scale T* 
in the e d /U vs Ae/U plane for V/t = 0.12, U/(2irt) = 1.0, 
and t' = t. The results in the range of -30 < log 10 (T*/T) < 
are painted in the colors varying from blue to red. T* 
becomes very small in the purple region, and the solid lines 
there denote the contours for log 10 (T*/T) = —30, —60, and 
—90. The dashed lines denote the phase boundary in the limit 
of T — >• 0, corresponding to Fig. [6] (b) . 



such that e apcx < Ep and e^ > Ep, with the Fermi energy 
Ep = 0.0. Therefore, the single electron enters mainly 
the apex site, and the other two dots are almost empty. 
This can also be confirmed through the charge distribu- 
tion plotted with the dashed line labelled "NlS05:even" 
in Fig. Q1JJ For Ae < 0.0, the occupation number 7V Q o for 
the apex site approaches to 1.0 while that for the even a\ 
orbital, N a i, almost vanishes. Therefore, the screening of 
the local moment is achieved through a super-exchange 
process by the conduction electrons which come to the 
apex site over the potential barrier at the other two dots, 
and thus T* decreases in this parameter region. 

Similarly, T* is suppressed also in the five-electron re- 
gion for Ae > 0.0, which corresponds to the area for 
£d/U < -1.2 at the top left of Fig.HH] The ground state 
in this parameter region has an even parity, and in the 
limit of T — s- the eigenvector is given by | $iven) in Eq. 
(IB13p . The average number of electrons N a0 and iV a i 
for this state coincide with those with respect to |$„dd) 
defined in Eq. (|B12I) . The results are shown in Fig. [T9l 
with the dot-dash line labelled "N4Sl:odd". We can see 
that a single hole with a spin 1/2 stays in the apex site 
for Ae > 0.0, and the other two dots are almost dou- 
bly occupied. Therefore, T* becomes very small in this 
case. This can also understood from the feature of the 
spectral function, shown in Fig. 1201 The Fermi level for 
the five electrons in this case is situated in the middle of 
the sharp peak at w ~ 1.2i in Fig. [20] (a). The width of 
this resonance corresponds to T* in the noninteracting 
limit, and the peak will become much narrower for finite 
Coulomb interaction U. There is another five-electron re- 
gion for A < 0.0, where the eigenvector is given by l^^d) 
in Eq. (|B14|) for the isolated TTQD. The local moment 
in this case stays at the b\ orbital, which is close to one of 
the leads, and thus the screening can be achieved by the 




(a) Ae/U (b) Ae/U 



FIG. 19. (Color online) Charge distribution in the isolated 
TTQD, N a ,i = ~Yli a ( na - ia ) > i s plotted as a function of Ae/U for 
(a) N a o the apex site, and for (b) the even ai orbital described 
in Appendix [B] The average is taken with respect to some of 
the eigenstates of Hdot = Udot + ^dot; an< l labelled by the 
occupation number, spin, and parity ("N4Sl:odd" denotes, 
Ntot = 4, S = 1, and an odd parity). The parameters are 
chosen to be T = 0, U/(2ivt) = 1.0, and t' = t. Note that 
Ae = e ape x — td, and the occupation number takes the form 
Mot = N a , + N a>1 + N b ,i, 



usual kinetic exchange mechanism of the 5 = 1/2 Kondo 
effect. 

The screening temperature becomes small also in the 
three region, and in the four electron region. In the three 
electron region for Ae > 0.0 the ground state has an odd 
parity, and the charge distribution for this state is plotted 
in Fig.m with the dotted line labelled "N3S05:odd" . We 
can see that at Ae ~ 1.0 the three electrons distribute 
almost homogeneously as N aQ ~ N al ~ 1.0, and thus 
Nf,i ~ 1.0. Then, as Ae increases, a single electron in 
the apex site moves towards the even-parity a\ orbital, 
and the occupation numbers approach to N a o ~ 0.0 and 
N a i — 2.0, keeping the occupation of the odd-parity or- 
bital almost unchanged N i ~ 1.0. Therefore, T* be- 
comes larger as Ae increases. 

The ground state in the other three electron region, 
for Ae < 0.0, has an even parity. The charge distribu- 
tion with respect to this state is shown with the solid 
line labelled "N3S05:even" in Fig. [Ill and in this case it 
is such that N a0 ~ 1.0, N ai ~ 1.5, and Nbi — 0.5 near 
Ae ~ 0.0~. Therefore, the apex site is singly occupied, 
and the moment protracted. The distribution varies as 
Ae decreases, and the fraction of the local moment tends 
to stay close to the leads, as iV a o ~ 2.0, N a i ~ 0.5, and 
N i ~ 0.5. We have also confirmed that T* is enhanced 
near Ae ~ 0.0 and —0.7 < ed/U < —0.4, along the sharp 
conductance valley caused by the SU(4) Kondo effect. 
The local moment in the three-electron region also re- 
flects the properties of the one-particle state. Similar 
to the t'/t > 1.0 case discussed in the previous subsec- 
tion, the orbital degeneracy is lifted as EP > E^ } + for 
Ae < 0.0. In this case the third electron enters the reso- 
nance state corresponding to the excited even-parity or- 
bital, appearing in the middle of Fig. (b). Since the 
spectral weight of this state is mainly at the apex site, the 
resonance width becomes narrow already in the noninter- 
acting case and it evolves into a sharp Kondo resonance 
for finite U. 
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FIG. 20. (Color online) Spectral functions in the noninter- 
acting case for (a) Ae = 0.5t and (b) Ae = —0.5t. The 
solid line is the even component — Im [G a o(u) + G a i(w)]/7r, 
and the dotted line is the odd component — ImG6o(<^)/T, de- 
fined in Appendix [B] The other parameters are chosen to be 
t'/t = 1.0, e d = 0.0, and T/t = 0.12. 



The S = 1 Kondo effect takes place in the four electron 
region at -0.4 < Ae/U < 0.8 in Fig. [HI We can see 
that T* varies significantly inside this region depending 
on whether Ae < 0.0 or Ae > 0.0, although the S = 1 
Nagaoka high-spin state which has an odd parity evolves 
continuously as Ae varies. The eigenvector in the limit 
of r — >• is given by |$ot d ), and the charge distribution 
with respect to this state is shown in Fig. [HO with the 
dot-dash line labelled "N4Sl:odd". We can see that N a0 
and N a i vary rapidly near Ae ~ 0.0, keeping the filling of 
b\ orbital unchanged N\>\ = 1.0. For Ae > 0.0, the local 
moment has a finite component in the apex site. This 
component of the moment moves to the a\ orbital near 
the leads for Ae < 0.0. The variation of T* inside the 
5 = 1 Kondo region reflects these changes in the charge 
and spin distributions. 



VII. SUMMARY 

We have studied the effects of distortions which break 
the full symmetry of an equilateral triangle of a TTQD 
connected to two non-interacting leads, over a wide range 
of the gate voltage ed- Two types of disorder have been 
considered, (i) an inter-dot tunneling matrix element t' 
t), and (ii) a level position e apcx e^) of the dot at 
the apex site. We have concentrated on the low energy 
behavior, restricting attention mainly to the regime with 
large Coulomb interaction U and small hybridization T 
as this leads to several different types of the Kondo effect. 

We find that the key variables for characterizing the 
low energy behavior are the total occupation number 
iVtot and the phase difference = (S c — 5 )(2/ir). The 
two phase shifts for the renormalized quasi-particles, 5 e 
and 5 , can be deduced theoretically from the low en- 
ergy NRG fixed point. The phase shifts may be deduced 
experimentally through the conductances and TVtot- Mea- 
surements of the AB oscillation in a magnetic field may 
also give a clue to determine the phase difference. 

In the parameter space for large U we find plateaus 
with the integer values of 0, and at each plateau the 
occupation number also approaches to an integer. These 



plateaus, therefore, can be classified with the two integer 
set (iVtot, O) [see Figs. [TT1 and [T4"]. The structure of these 
plateaus of determines the precise feature of the Kondo 
ridges and valleys of the conductance [see Figs. 151 and[T3"]. 

Different Kondo effects occur in different regimes. The 
SU(4) Kondo effect takes place for weak distortions, 
along the contour line for = 2.0 which runs in the 
region of iVtot — 3.0 in the parameter space. This con- 
tour transverses the middle of a steep cliff of 0, standing 
between the plateau for ~ 1.0 and that for ~ 3.0. 
It can be observed as a sharp conductance valley be- 
tween the Kondo ridges on both sides, and the slope of 
the cliff determines the width of the conductance valley. 
The SU(4) Kondo behavior is sensitive to the perturba- 
tions which lower the symmetry of the equilateral trian- 
gle. This is caused by the fact that the SU(4) symmetry 
relies crucially on the orbital degeneracy. Furthermore, 
the spin and charge distributions inside the TTQD vary 
near the SU(4) symmetric point, and it affects signifi- 
cantly the Kondo screening. 

The S = 1 Kondo effect, taking place at the plateau of 
for (iVtotjO) — (4.0, 2.0), is robust against the break- 
ing of the symmetry of the equilateral triangle. This is 
mainly due to a size effect: there is a finite energy separa- 
tion between the Nagaoka high-spin state and the excited 
local singlet state in the isolated TTQD cluster. For large 
distortions a singlet-triplet transition takes place. It be- 
comes a crossover between a Kondo and non-Kondo sin- 
glet state for finite T, and the series conductance has a 
peak of the height of 2e 2 /h in the transient regions. The 
width of the peak is determined by the slope of the cliff 
of 0, which appears at the crossover region. 

Apart from the phase shifts which determine the con- 
ductance and occupation of the TTQD, another impor- 
tant renormalized parameter characterizing the low en- 
ergy behavior is the characteristic energy scale T*. For 
T < T* the low-energy properties can be described by the 
local Fermi-liquid theory. In the cases where the quan- 
tum dots have a local moment, T* can be regarded as 
the Kondo temperature. We have estimated T* from the 
region where the NRG levels crossover to the low energy 
fixed point. The results for T* reflect the distribution of 
the charge and spin in the TTQD [Figs. [151 and ITS]. 

Specifically, T* tends to be small in the case where a 
partial moment remains in the apex site, which has no 
direct coupling to the leads. The screening of such a 
partial moment becomes sensitive to the charge and spin 
on the other two dots because the conduction electrons 
tunneling from the leads have to pass through cither of 
the two dots to get the apex site. In some regions of the 
parameter space, we find that the tunneling of the con- 
duction electron is suppressed at these two dots, in a way 
analogous to a super-exchange process caused by a po- 
tential barrier between the local moment and leads. The 
characteristic temperature T* can be raised, however, by 
making the coupling to the leads T stronger. Note that 
T* depends on T not only through the prcfactor, but also 
through the higher order contributions of the hybridiza- 
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tion, which cause an exponential dependence of T* on T 
and other parameters. Specifically, T* may become large 
for the TTQD with a small charging energy U. Our re- 
sults provide an overview of how characteristic energy 
scale varies in the different the regions in the parameter 
space. 

A general point worthy of note is that the two types 
of the distortions show a clear contrast in the form of 
the charge distribution for some regions of the parame- 
ter space [see Figs. [TBI and [19] . The diagonal distortion 
(fapcx 7^ e<i) affects directly the potential of the apex site, 
so that the charge distribution is more sensitive to Ae 
than to the off-diagonal one (f ^ t), and this difference 
affects the characteristic energy scale T* significantly in 
some regions of the parameter space. 
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FIG. 21. Even and odd orbitals: the onsite potential of each 
orbital is given by e a0 = e apC x, e a i =£<( — *', and eu — e d + t'. 



at the Fermi level u = 



26.28 



(A3) 



-r L ImG+(0)-r fl ImG+ DWD (0) . (A4) 



g s = ^-4T B r L \G+ Dl (0)\ 2 
2e 2 r 

Note that the contributions from the vertex correction do 
not appear here due to the property that the imaginary 
part of the self-energy vanishes Im (0) = at T = 
and u> = 04£ Furthermore, for the symmetric coupling 
Ti = Tr (= r), the Green's functions can be expressed 
in the forms, 



G+(0) = G 



N D N, 



,(0) 



1 

2T 



1 



1 



k c + i 
1 

k c + i 



K + I 
1 

Ko + i 



(A5) 
(A6) 



Here, n c = — cot S e and k q = — cot S a include all the 
many-body corrections, through the real part of the self- 



energy ReEy(0)i=a Equations 
1. 



follow from Eqs. 



Appendix A: Phase shifts S c and S 



Appendix B: Even and odd orbitals 



The phase shifts for interacting electrons can be de- 
fined, using the Green's function 

Gij(uj n ) = - fdr (r r ^(r)4(0)) e^ T . (Al) 



Here, = 1/T, dUr) = e^d^e'^ , and (O) = 
Tr [e-P n O] /Tre~P n . The retarded Green's function is 
given by Gfj(uj) = Gij{u) + i + ) via the analytic contin- 
uation, and the self energy Eij(z) due to the interaction 
H^ ot can be described by the Dyson equation 

N D N D 

G l0 (z) = G[f(z) + E E G u'(z) Ki'Az) G f3 {z) . 

i'=l j' = l 

(A2) 

Here, the number of the dots is Nd = 3 for the TTQD, 
and Gf){z) is the non-interacting Green's function cor- 
responding to the free Hamiltonian Ho = W° ot + Hmix + 

At zero temperature T = 0, the series g s and parallel 
g p conductances are determined by the Green's functions 



The eigenstates of the Hamiltonian T~L defined in Eq. 
(P) can be classified according to the parity in the case 
that the system has an inversion symmetry, using the 

d 



even-odd basis defined by a 0a = u 2a , 



'he 



V2 



he 



(Bl) 



The labels and 1 for the even-odd basis are assigned 
in the way that is shown in Fig. 1211 The odd parity b\ 
orbital corresponds to the eigenstate for Eo , defined in 
Eq. (fT5|) . for the noninteracting TTQD cluster. Similarly 

the eigenstate for is given by a linear combination 
of the even ao and a\ orbitals. 

|$«) = (u+al + u.al) |0), (B2) 
|$«) = (-u_aL + u + a\ a ) 10) • (B3) 

Here, |0) is a vacuum. The coefficients are normalized 
such that u + + u 



1. 



1 ± 



Ae + t' 



(B4) 
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and u 2 , {u 2 _) corresponds to the spectral weight for the 

a (m) component in the excited state l^i)- When the 
TTQD is coupled to the leads, these states become the 
resonance levels, which can be described by the Green's 
functions for the do, <Zi, and b\ orbitals 



1 



2t 2 



capex u-Ea+t'+iT 
1 



UJ — €d + t' + iT 

1 

uj — e c i — V + iT 



(B5) 
(B6) 
(B7) 



Specifically for Ae + 1' > the coefficients take the value 
of u\ > 0.5 and u 2 _ < 0.5. In this case the even excited 

state becomes a sharp resonance peak, the spec- 

tral weight of which is mainly on the apex site, and the 
spectral function near u) ~ E^ 1 } takes the form 



-ImGi^w) * u 



uj- E, 



ImG^) ~ u 2 _ 



(i) 
u 2 _Y 



u 2 _T) 



uj-E. 



(i) 

c + 



, (B8) 



(B9) 



u 2 _r) 



The Green's functions near the lower level uj ~ E c _ can 
also be written in similar forms, just by replacing + (— ) 
in the suffix by - (+) in Eqs. JBSJ} and (|H5)l . 

The interaction Hamiltonian defined in Eq. (J3)) can be 
expressed, in terms of these even-odd orbitals, in the form 

%dot = U n a ,of n aM + — (n a s t n a ,n + «MT n b,n) 

+ j (^2 na l nb ' 1 ~ 2 '^ a ' 1 ' ^ b < 1 

+ \ ( a n a n b u b n + b n b u a n a it) ■ ( B1 °) 

Here, S a ,i = a \a&<?<?' a i<7' A ° tnc Pauli matrices, 

n a ,i<j = a\ a a icn and n 0)i = J^a n a,ia- The operators n bM 
and 5t,,i for the odd-parity orbital b\ a are defined in the 
same way. 

The Hilbert space for the isolated TTQD cluster, de- 
scribed by "H dot = %°ot +^doti can be constructed from 
the three orbitals. For instance, the Nagaoka state for 
H dot has an odd parity. It has one electron in the b\ 
orbital, and the eigenvector takes the form 



l^dd)= «o|I)+ai|H>, (Bll) 
'I^ a it a ot a o4,|0)> I- 1 - 1 / 



|I) = 6+ aLoLat |0), |II) = 6LoLoI,oL|0) . (B12) 



Here, ao and ot\ are the coefficients. The eigenvectors for 
five electrons can be expressed in the form 



|$(5) 
I evei 

(5) 



"Ul^odd/ 



(B13) 

^lt^lt^H^Ot^OJ.!"/- ( B14 ) 

Therefore, the odd-parity b\ orbital is fully occupied for 



1$ 



odd/ 



ftLoLalaLaJjlO). 



I* 



(5) 



while the even ao and a\ orbitals are fully occu- 



pied for |$ odd ). The distribution of the charge and spin 
in the three orbitals affects significantly on the way the 
screening by the conduction electrons is carried out for 
finite r, when the leads are connected to the TTQD. 



Appendix C: NRG approach 



We provide an explicit form of the discretized Hamil- 
tonian Hn of the NRG in this appendix. The non- 
interacting leads are transformed into the tight-biding 
chains in the NRG approach, through the logarithmic 
discretization with the parameter A. Then, a sequence 
of the Hamiltonian Hn with a finite size is introduced in 
the formic 



Hn =A (N - 1)/ Unl ot +H^ t + H mU + H^i), (CI) 



'0, La 



H mix = «E(/<UA* + d{j { 

a 

+ (foMa d N c ,a + d N a ,Jo 



Ra I 1 



(C2) 



w 1 + VA 



N-l 



load 



E EE^«a 

v=L,R a n=0 



-n/2 



x ( fn+l,vo /n,i/<T + fn,va fn+l,L>a ) ■ (C3) 



Here, D is the half-width of the conduction band, and 
the other parameters are defined by 



2DTA A 1 1 + 1/A 

~^T~ ' Aa ~ 2 i~T/X 

1 - 1/A™ +1 



y/i - i/A 2 "+yi - 1/A 



2n+3 



log A, (C4) 
(C5) 
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